Decay of nuclear hyperpolarization in silicon microparticles 
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We investigate the low-field relaxation of nuclear hyperpolarization in undoped and highly doped 
silicon microparticles at room temperature following removal from high field. For nominally undoped 
particles, two relaxation time scales are identified for ambient fields above 0.2 mT. The slower, Ti^s, 
is roughly independent of ambient field; the faster, Ti,f , decreases with increasing ambient field. A 
model in which nuclear spin relaxation occurs at the particle surface via a two-electron mechanism 
is shown to be in good agreement with the experimental data, particularly the field-independence 
of Ti^s. For boron-doped particles, a single relaxation time scale is observed. This suggests that for 
doped particles, mobile carriers and bulk ionized acceptor sites, rather than paramagnetic surface 
states, are the dominant relaxation mechanisms. Relaxation times for the undoped particles are not 
affected by tumbling in a liquid solution. 



I. INTRODUCTION 

Silicon has long been a staple of the microelectron- 
ics industry and so has been the subject of intense ma- 
terials research. In recent years, nanoscale silicon has 
attracted attention due to its unique electronic and op- 
tical properties, and as a potential agent for medical 
imaging and drug delivery.'^"Sl Spin phenomena in sili- 
con has been investigated using nuclear magnetic res- 
onance (NMR) for over a half-centuryj^^Ki^ while mate- 
rial properties and fabrication has continually developed. 
For instance, the nuclear spin-lattice relaxation tim e, Ti, 
provides information about dopants and impurities .Hm^l 
Notably, the low natural abundance (4.7%) of spin- 
1/2 ^^Si nuclei in a lattice of zero spin nuclei leads to 
Ti of many hours® ^'^ and coherence times up to tens of 
secondd^ in undoped bulk crystalline silicon. 

These remarkable NMR properties have stimulated in- 
terest in silicon as a platform for solid-state quantum 
comput ing liSHiZI ^nd as a long-lived imaging agent for hy- 
perpolarized magnetic resonance imaging (MRI).'S'In par- 
ticular, there have been renewed efforts to understand dy- 
namic nuclear polarization (DNP), a process where sat- 
uration or pumping of paramagnetic impurity states by 
microwave fields or optical illumination can lead to nu- 
clear spin polarizations orders of magnitude larger than 
thermal equilibrium polarizations.^® Enhancing polar- 
ization by these methods is commonly referred to as hy- 
perpolarization. 

Application of silicon hyperpolarized imaging 

agent for MRI requires an understanding nuclear spin 
relaxation over a broad range of applied fields, from 
high fields where polarization is induced and imaged, to 
low fields where the agent is transferred from the po- 
larizer and administered. Previous detailed studies have 
focussed on low-temperatures, investigating the depen- 
dence of Ti on applied magnetic field, dopin^^MUl and 
strain. lafil It was found that Ti increased with applied 
magnetic field, which was explained in terms of the field 
dependence of electron- nuclear dipole-dipole interaction. 



The nuclear Ti was also shown t o scal e inversely with the 
concentration of mobile carriers .'i^'^ In silicon nanopar- 
ticles, the observed increase in nuclear Ti with increasing 
particle size was attributed to diffusion-mediated relax- 
ation via defects at the particle surface.!^ Recent mea- 
surements have shown that the relaxation and coherence 
times of "^^P donor-bound electron spin are dramatically 
reduced near the silicon surface due to Pb defects, in 
comparison to known values in the bulk.'^ Studies of 
electrically^ and optically— detected hyperpolarization 
of ^^P donor nuclear spins near the silicon surface, how- 
ever, do not consider effects of the surface on nuclear spin 
relaxation. 

In this paper, we investigate the decay of hyperpolar- 
ization of ^^Si nuclei in silicon microparticles at room 
temperature and low ambient magnetic fields after re- 
moving the samples from a high-field polarizing environ- 
ment. We find the decay is bi-exponential, with a slow 
time scale that is independent of ambient field, and a 
fast time scale that shows only a modest decrease with 
increasing ambient field. We develop a model of nuclear 
spin relaxation in silicon particles that takes into account 
the heterogeneous makeup of the sample, with direct nu- 
clear relaxation occurring only near the surface, while 
nuclei in the particle core are relaxed indirectly by spin 
diffusion. This model extends a previous spin diffusion 
model, ^® which predicts Ti increasing with the square of 
the particle diameter. We find that the weak magnetic 
field dependence observed experimentally is inconsistent 
with a simple extended spin-diffusion model based on re- 
laxation on individual bound electrons. However, by gen- 
eralizing the model further to include nuclear spin relax- 
ation mediated by pairs of dipolar-coupled electrons, the 
essentially field-independent relaxation is recovered by 
the model. "^^ 33 35 experimental methods presented 
here, as well as the use of magnetic field dependence of 
nuclear Ti as a probe for testing the relative importance 
of these different model processes, have precedent in a 
similar investigation performed on diamond and its ^^C 
nuclear spin relaxation .I^SHMl 
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II. SAMPLES AND MEASUREMENTS 

Microparticles made from nominally undoped and 
boron-doped ball- milled Si wafers were investigated. As 
shown previousljS, undoped Si microparticles have very 
long high-field relaxation times (Ti ~ 10^ s), relevant for 
hyperpolarized MRI. Undoped samples were produced 
by ball-milling high-resistivity (> lOkSl-cm) float-zone 
Si wafers (Silicon Quest International) followed by size 
separation by centrifugal sedimentation, yielding a mean 
particle diameter of 5 /im.^ Boron-doped Si micropar- 
ticles with doping density ~ 5 x 10^^ cm^'^ have much 
faster nuclear spin relaxation, Ti 10^ s at high mag- 
netic field. The boron-doped microparticles were pro- 
duced by ball-milling Czochralski-grown Si wafers (Vir- 
ginia Semiconductor, resistivity 0.01-0.02 Jl-cm) followed 
by size separation, yielding a mean particle diameter on 
the order of 10 ^m. Except where noted, samples con- 
sisted of a 0.1 mL teflon tube filled with particles. All 
NMR measurements were carried out at Bq = 2.9 T using 
a custom-built probe and spectrometer. 

Polarization at 2.9 T following a saturation sequence 
of sixteen tt/2 pulses was measured function of 

time, ipoh using a Carr-Purcell-Meiboom-Gill (CPMG) 
sequenccPHm^ - [t/2 - (tt)^ - t/2 - echo]" with 

t — 1 ms and n — 400. After each data point, the sam- 
ple was re-saturated and the measurement repeated. As 
seen in Figs. 1(a) and 4(a), the build-up of polarization 
is well described by a single exponential function of tpoi 
for both undoped and doped samples. 

Depolarization at low ambient fields was measured af- 
ter first polarizing at 2.9 T for 8 h (- 3Ti at 2.9 T), 
then raising the probe out of the magnet bore to a posi- 
tion where ambient (or holding) fields of 0.2 mT, 6 mT, 
130 mT, and 300 mT had been previously calibrated us- 
ing a Lakeshore 460 gaussmeter. In addition, a nomi- 
nal zero-field measurement used a commercial zero-gauss 
chamber (Lakeshore 4060) with ambient field below 1 /iT. 
Following depolarization in ambient field, the sample was 
returned to field center (2.9 T) and the remaining nuclear 
polarization was measured using the CPMG sequence de- 
scribed above. The transit time between positions was 
~ 5 s, fast compared to Ti but slow compared to nuclear 
Larmor times. The sample was fully repolarized follow- 
ing each data point. 

In addition to NMR measurements, room temperature 
electron spin resonance (ESR) measurements were car- 
ried out to characterize electronic defects and carriers of 
the two sample types. The ESR spectrum of the undoped 
sample showed a single peak at _B = 324 mT measured 
at ESR frequency / = 9.099 GHz, corresponding to a 
g-factor centered at 2.006, with a linewidth of 0.47 mT. 
The ESR spectrum of the boron-doped sample showed 
a single peak at B = 336 mT measured at ESR fre- 
quency / = 9.444 GHz, again corresponding to a g-factor 
of 2.006. A g-factor of 2.006 is consistent with reported 
values for Pb defects at the Si-Si02 interface, and is not 
consistent with g-factors of other common dopants or ox- 
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Figure 1 : Relaxation of nuclear spin polarization of undoped 
silicon particles, (a) Recovery of spin polarization following 
saturation at Bo = 2.9 T. (b) Decay of nuclear magnetization 
at -Bdcp = 0.2 mT (circles), 6 mT (open squares), and 300 mT 
(triangles), after polarization at 2.9 T, and same data for par- 
ticles suspended in methanol at 6 mT (filled squares). Curves 
show best fits to Eq. 1. No signal was observed after 1 second 
for data taken in the zero-field chamber (< 1/iT). 

ide states.ESHHl 



A. Undoped Microparticles: Measurements 

Polarization and depolarization data for the un- 
doped microparticles are shown in Fig. 1. Po- 
larization is well described by a simple expo- 
nential, P ^ fb,Bo (l - "^i*^^ best-fit 
value Ti = 8200 ± 600 s, consistent with previ- 
ous measurements.'^ Equilibrium polarization is small, 
Po.So = tanh[{h-fnBo)/{2kBT)] = 2 x 10-^. Here, ke 
is Boltzmann constant, T = 300 K is room temperature, 
and 7„ = 5.31 x 10^ T^^ is the nuclear gyromagnetic 
ratio for ^^Si. 

In contrast to the build-up of polarization at 2.9 T, 
depolarization at low ambient fields, -Bdcp- decays with 
two distinct time scales, which we characterize using a 
bi-exponential form, 

P = Po.B,^, + {Po.Ba - Po,B,J (1) 

X (^{1 ~ Q;)e^*'^'=p/^i'f -I- ae"*<i''p/'^i-=^ , 

where Ti j and Ti s are fast and slow relaxation times, and 
a is the fraction of spins whose polarization decays slowly. 
Best-fit values of a, Ti^f and Ti.g are shown in Table |T] 
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Bdcp 


a 


Ti,t (s) 


Ti.s (s) 


300 mT 


0.765 ± 0.027 


531 ± 179 


13850 ± 1030 


6 mT 


0.593 ±0.050 


903 ± 232 


15310 ±2070 


0.2 mT 


0.591 ±0.120 


1124 ± 544 


14910 ± 5360 


< 1 fj,T 




< 1 





Table I: Exponential weights and relaxation times fitting to 
experimental data. Ti^f and Ti,s are fast and slow relaxation 
times with which two additive components of the spin polar- 
ization decay, a is the relative amplitude of the component 
which decays at Ti,s. 



Note that the slow relaxation time depends very weakly 
on Bdcp, with Ti^s ^ 1-4- 1.5 x 10'^ s in all cases. In con- 
trast, fast relaxation becomes somewhat faster (shorter 
Ti f) with increasing -Bdcpi while the fraction of slow re- 
laxers increases. 

Inside the zero-gauss chamber the residual ambient 
magnetic field is specified to be below the few-/iT scale, 
much weaker than the typical nuclear dipolar field, 
-Sn.dd /^oMn/(47ra'^) = 0.08 mT, where a ~ 4 A is the 
mean separation between randomly distributed ^^Si nu- 
clei. At such extremely low fields, spin transitions will oc- 
cur that are forbidden when B > i3„ by the conservation 
of nuclear Zeeman energy, and nuclear spin relaxation is 
as fast as decoherence.^E As expected, we find that in 
the zero-gauss chamber nuclear polarization decays very 
quickly, with no visible signal for i^ep > 1 s. 

To investigate the effect of microparticle tumbling on 
depolarization, a fluid mixture of one part Si micropar- 
ticles and four parts methanol by weight was compared 
to a sample of packed dry powder. Depolarization times 
at 6 mT ambient field showed no significant difference in 
relaxation times. 
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B. Undoped Microparticles: Model 

We interpret NMR and ESR data for the undoped mi- 
croparticles within a model shown schematically in Fig. 2, 
comprising nuclear dipolar diffusion in a core region and 
relaxation via paramagnetic sites in a shell region near 
the Si/Si02 interface. In undoped silicon, direct cou- 
pling of ^^Si nuclear spins to phonons is weak^i*^, and 
dipolar coupling to paramagnetic impurities and defects 
typically dominates nuclear spin relaxation. Within our 
model, the bi-exponential form of Eq. 1 refiects the ex- 
istence of two populations of nuclear spins within each 
microparticle: a fraction a, located in the core of the 
particle, has a long spin relaxation time Ti^g mediated 
by nuclear spin diffusion; the remaining fraction (1 — a), 
located within a shell near the Si/Si02 interface, has a 
short relaxation time Ti f mediated by electron- nuclear 
spin interaction associated with paramagnetic centers at 
the Si-Si02 interface. 

For i?dcp > 5n,dd! the nuclear dipolar spin diffusion 
rate is independent of magnetic field,'^ and is well de- 
scribed by a diffusion constant D = Wa^ ~ a^/(50r2), 
where W is the probability of a flip-flop transition be- 
tween nuclei due to dipole-dipole interaction, and T2 is 
the nuclear decoherence time.'^'^ This is consistent with 
our observation that Ti.g is roughly independent of -Bdcp 
across a broad range of values. 

The thickness of the shell region is set by the nuclear 
spin diffusion radius, /3 — {C/DY^'^, where C is a con- 
stant describing the nature of the dipolar interaction oc- 
curring between the nuclei and electrons. Nuclei situated 
at a distance r from the electrons at the particle surface 
may relax through a dipolar interaction with the elec- 
trons at a rate cx C/r^^ 

We consider two physical mechanisms of relaxation in 
shell region, one in which each nuclear spin interacts with 
a single electron spin, and the other in which each nuclear 
spin interacts with a pair of electron spins. The two- 
electron model captures a crucial physical mechanism by 
allowing energy matching between electron-spin-pair flip- 
flops and nuclear spins flips in a magnetic field, leading 
to efficient nuclear relaxation with weak dependence on 
magnetic field. 

Modeling nuclear spin relaxation as arising from inter- 
action with individual paramagetic centers is appropri- 
ate for sparsely distributed impurities. This regime has 
been studied experi menta lly, for instance, in diamagnetic 
crystals such as LiJ^^IIill a,nd CaF2.'^ In this case, the 
orientation-averaged dipole interaction strength is given 

by 



Figure 2: Model showing the spin reservoirs and relaxation 
pathways in undoped silicon particles. Nuclei in the core re- 
gion of the particle relax by transferring their magnetization 
to the nuclei in the shell region, nearer to the surface, by 
spin diffusion. The nuclei near the surface can relax quickly 
due to strong dipolar coupling to paramagnetic defects at the 
Si/Si02 interface. 



C(le) = 



10 TieBl 



(2) 



dcp 



where 7e is the gyromagnetic ratio of the electron and Tie 
the electron spin-lattice relaxation time.'22l This model 

yields a magnetic field dependence of diffusion-limited 

1/2 

relaxation (i? ^ /?) at least as strong as Ti s oc B^^^ 
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because the diffusion radius /? and relaxation rate cx 
C depend explicitl y on field, even when D is field 
independent This prediction is incompatible with 
our experimental results, which show a much weaker field 
dependence. 

In light of this inconsistency, we consider a model of 
nuclear relaxation at the particle surface that includes 
three-spin processes involving pairs of interacting elec- 
trons, which gives 



C(2e) - 



3 



2^,2 



^dcp7nT'2e 



n2 ^27^2 



dw, 



(3) 

where T2e is the electron spin-spin coupling and 
g{oj) is the normalized electron absorption lineshape 
function .'^flE^ This model has been previously applied to 
systems with more concentrated par amagn etic impuri- 
ties, including La2Mg3(N03)i2 • 24H2Cl^^an d mor e re- 
cently ^^C nuclear spin relaxation in diamond .I^SEI] This 
model accounts for flip-flop transitions between nearby 
electron pairs, occurring on a time scale Tie, which 

provide the fluctuating magnetic fleld that can flip nu- 
clear spins. That is, when the dipolar coupling of elec- 
tron pairs matches the nuclear Zeeman energy, a three- 
spin interaction can occur that exchanges the spins of the 
two electrons while flipping a nuclear spin.SslThis process 
depends on the density of transitions between electronic 
dipolar energy states that match the nuclear Zeeman en- 
ergy. For a Lorentzian electron lineshape, the integral in 
(3), which describes the probability of flnding two elec- 
trons within the ESR line differing in frequency by a;„, 
can be replaced by 2/(4 + a;^T|g). For low B^cp, such 
that 7„i3dcp < leBcdd ^ T^ l, where i3o,dd is the elec- 
tronic dipole field felt by a typical nucleus, the density of 
nuclear-spin-fiip transitions will be independent of -Bdop, 
hence the rate of direct nuclear spin relaxation by this 
process will also be independent of Bdep-'^^'^ 

We infer an upper bound for i?c.dd from the measured 
ESR linewidth of 0.47 mT, which gives 7^^7e-Bo,dd ^ 
1.7 T for the scale of i?dop below which nuclear Ti 
should be roughly field independent within this three- 
spin model. The low end of the field-independent range 
is set by the nuclear dipolar field i?n,dd ~ 0.08 niT. 
This range is consistent with the experimental observa- 
tion that Ti s is roughly field independent from 0.2 mT to 
300 mT. Weak field dependences of Ti f and a are likely 
due to the weak field dependence of Tie. 

Very close to each paramagnetic impurity, nuclear spin 
diffusion is suppressed, creating a barrier to diffusion of 
radius b = a(fi.7gi?dep/7n2fcBTV/'*, due to gradients in 
the nuclear Larmor frequency.1221 This effect is included 
in the simulations (described below), but because (3 ^ b 
it only weakly affects overall relaxation times. 

We have simulated nuclear spin relaxation as a func- 
tion of ambient field for a spherical silicon micropar- 
ticle with native ^^Si concentration and paramagnetic 
sites at the surface. For both one and two electron pro- 
cesses, the average polarization P(idGp) is approximately 
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Figure 3: Experimental slow (black circle) and fast (black 
triangle) relaxation times, Ti^s and Ti^f, from Table [l] for 
bi-exponential nuclear spin relaxation, along with simulation 
results for a model that includes either one-electron or two- 
electron processes. 



bi-exponential, similar to the experimental data. Fits to 
the relaxation curves gave Ti^s and Ti^f values shown in 
Fig. [3j The particle diameter was taken to be 700 nm, 
smaller than the actual mean particle size. This com- 
pensated the spherical assumption, which gave too small 
of a surface where defects reside. Other parameters in 
the model were nuclear T2 of 10 ms (and correspond- 
ing nuclear spin diffusion constant 0.3 nm^/s), electron 
Tie of 200 ns, and electron of 25 ns. Figure |3] shows 
the striking contrast in field dependence between the two 
models, and that the three-spin process compares rela- 
tively well with experiment. 



C. Results for Boron-doped Microparticles 

Nuclear polarization and depolarization of boron- 
doped Si microparticles are shown in Fig. 4. Both P(ipoi) 
and P(idep) are well fit by single-exponential relaxation 
curves with time constants Ti that increases with increas- 
ing ambient field. The reduction of Ti by 1-2 orders of 
magnitude in the doped case reflects the dominant con- 
tribution of nuclear relaxation from mobile carriers,!^ 
as well as from ionized acceptors distributed uniformly 
through the particle instead of just at the surface (Fig. 5), 
bypassing the relatively slow nuclear spin diffusion pro- 
cess. For the specifled dopant density, the average dis- 
tance between dopants is ^ 2.5 nm and the Bohr radius 
^ 1.3 nrrPSl, both an order of magnitude less than the 
diffusion radius for both single electron and two-electron 
processes across the range of ambient fields under study. 
The contribution of mobile carriers can be expressed as 
an additional relaxation rate which is independent of 

oc NqA^ where iVo is concen- 
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Figure 4: Nuclear spin lattice relaxation in boron-doped sil- 
icon microparticles. (a) Recovery of spin polarization after 
saturation at _Bo = 2.9 T. (b) Decay of nuclear magnetization 
at -Bdcp = 6 mT (open square), and f30 mT (triangle), af- 
ter polarization at 2.9 T. Error bars in (a) are displayed but 
smaller than the data points shown. 



< — > paramagnetic centers 



< — lionized B acceptors 



► mobile liolesi 



direct relaxation 



Si nuclei in bull< 




Figure 5: Rate model showing dominant nuclear relaxation 
pathways of direct relaxation on mobile holes and short dif- 
fusion to ionized boron acceptors throughout the bulk of the 
particle. Long distance spin diffusion to surface impurities is 
relatively slow and direct coupling of the nuclei to the lattice 
is weak. 



III. CONCLUSION 

We have investigated the relaxation of polarization of 
natural-abundance ^^Si nuclei in undoped and doped sil- 
icon microparticle samples as a function of ambient mag- 
netic fields at room temperature. In undoped micron- 
scale particles, the decay of nuclear polarization from a 
hyperpolarized value was measured from microtesla to 
0.3 T. For fields stronger than the nuclear dipolar field, 
a bi-exponential decay of nuclear magnetization was ob- 
served, with a fast component of ~ 10^ s that depends 
weakly on field and a slow component of ^ 10^ s that 
is roughly field independent. The relative amplitude of 
the slowly decaying component increases with the mag- 
netic field, but only slightly. Bi-exponential relaxation 
suggests the presence of two nuclear spin baths distin- 
guished by their proximity to paramagnetic impurities at 
the particle surface. The timescales of nuclear spin relax- 
ation are largely independent of ambient field i?dop from 
0.2 to 300 mT. This independence from i?dcp is quantita- 
tively consistent with a model of nuclear spin relaxation 
dominated by a three-spin mechanism in which flip-flop 
transitions of two electrons at the nuclear Larmor fre- 
quency flip a nuclear spin. In highly doped silicon, sim- 
ple exponential relaxation with a faster Ti of order a few 
hundred seconds is observed. 

These results can be compared with similar experi- 
ments performed in diamond, another material with di- 
lute (1.1 %) spin-1/2 nuclei among a majority of zero-spin 
nuclei in an fee lattice-I^^El Measurement of the nuclear 
Ti in natural diamond as a function of magnetic field in- 
dicate a contribution of three-spin processes to ^■^C spin 
relaxation. '^'^ Other investigations of nuclear spin relax- 
ation in ^^C-enriched diamond^ and in nanocrystalline 
diamond"^^ have shown that the relaxation of nuclear 
magnetization was not well described by a simple ex- 
ponential approach to equilibrium. The NMR studies in 
nanocrystalline diamond revealed that the particles con- 
sist of a crystalline core and a surrounding shell. All 
these observations in diamond are similar to our results 
in silicon microparticles. 

Finally, we comment that multi-hour and field- 
independent nuclear relaxation times for a substantial 
fraction of the nuclear spin population, including for 
particles in a fiuid suspension, is important for for the 
use of silicon microparticles hyperpolarized MRI 

imaging agent. 



tration of mobile carriers, and A is the strength of the 
hyperfine coupling. The contribution to nuclear relax- 
ation by ionized acceptors as well as surface paramag- 
netic defects at the surface are responsible for the mild 
magnetic field dependence of Ti seen in Fig. 4(b). We 
were unable to observe an ESR signal from the mobile 
carriers, presumably due to their short relaxation times 
leading to broad broad ESR lines. 
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